Optical conductivity in the t- J-Holstein Model 
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Using recently developed numerical method we compute charge stiffness and optical conductivity of the t-J 
model coupled to optical phonons. Coherent hole motion is most strongly influenced by the electron-phonon 
coupling within the physically relevant regime of the exchange interaction. We find unusual non-monotonous 
dependence of the charge stiffness as a function of the exchange coupling near the crossover to the strong 
electron-phonon coupling regime. Optical conductivity in this regime shows a two-peak structure. The low- 
frequency peak represents local magnetic excitation, attached to the hole, while the higher-frequency peak 
corresponds to the mid infrared band that originates from coupling to spin-wave excitations, broadened and 
renormalized by phonon excitations. We observe no separate peak at or slightly above the phonon frequency. 
This finding suggests that the two peak structure seen in recent optical measurements is due to magnetic excita- 
tions coupled to lattice degrees of freedom via doped charge carriers. 



PACS numbers: 71.10.Pm,71.27.+a,78.67.-n,71.38.-k 
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I. INTRODUCTION 

Despite many years of intensive research of transport prop- 
erties of a hole doped in an antiferromagnetic background the 
proper description of this system remains a challenging the- 
oretical problem. The transport of a doped hole leaves in 
its wake locally distorted, slowly relaxing spin background, 
leading to the formation of a dressed quasiparticle with an 
enhanced effective mass and renormalized charge stiffness - 
a measure of a coherent, free particle like transport. Addition 
of lattice degrees of freedom to this already elaborate problem 
reflects the current scientific interest in the field of correlated 
electron systems. 

Long after the pioneering work 1 , the enhanced interest in 
correlated models, coupled to lattice degrees of freedom is pri- 
marily fuelled by experimental evidence given in part by an- 
gular resolved photoemission data demonstrating that strong 
electron-phonon (EP) interaction plays an important role in 
low-energy physics of high-T c material s 2 ' 3 ' 4 ' 5 . Moreover, re- 
cent estimates of transport based on the pure the t — J model 6 
yield substantially smaller resistivity in comparison to exper- 
iments in the low-doping regime of La2_£ C Sr a; Cu04, that can 
be explained as a lack of additional lattice degrees of freedom. 

Recent numerical methods investigating optical conductiv- 
ity (OC) in correlated electron systems and systems, where 
electrons are coupled to bosonic degrees of freedom, have 
been focused on the generalized t — J model 7,8 , Holstein and 
generalized electron-boson model 910 , while investigations of 
the t — J-Holstein model have been until recently limited 
to small clusters with 10 sites 11 . Diagrammatic Quantum 
Monte Carlo (DMC) method has been applied to resolve OC 
of the Frohlich polaron 12 and recently also the t — J-Holstein 
model 13 . In the later work authors report on a two-peak struc- 
ture in the optical response where the Iow-lj peak is due to 
polaronic effects while the peak at higher oj is due to magnetic 
excitations, renormalized by lattice degrees of freedom. DMC 
method does not reproduce the well known string states seen 



in the spectral function 14 neither the signature of local mag- 
netic excitations in optical conductivity since authors use the 
self consistent Born approximation (SCBA) without magnon- 
magnon vertex corrections for treatig spin degrees of freedom. 
Their result seems to contradict calculations based on the dy- 
namical mean-field theory (DMFT) where a two peak struc- 
ture is seen in OC only at small J/t < 0.3^. In this work 
authors also show that the low-a; peak is of the magnetic ori- 
gin while the higher-w peak represents the broad polaronic 
band. 



II. MODEL AND NUMERICAL METHOD 

The main goal of this work is to investigate in depth optical 
properties of the t- J-Holstein model for the case of a single 
hole in the antiferromagnetic background. We first define the 
t- J-Holstein model on a square lattice 

h = -t ]T c, s^ ■ ./^S,S.. 

(*>j}>« 

+ g^2(l-rii)(af + a,i) -{-vo^ afcii, (1) 

i i 

where £j )S = Cj )S (l — n,,_ s ) is a fermion operator, projected 
onto a space of no double occupancy, t represents nearest- 
neighbor overlap integral, the sum runs over pairs of 
nearest neighbors, a% are phonon annihilation operators and 
m = 53 n i-s- The third term represents EP coupling g = 
\/8\uJot, where A is the dimensionless EP coupling constant, 
and the last term represents the energy of Einstein phonons 

U>Q. 

We use recently developed method based on the ex- 
act diagonalization within the limited functional space 
(EDLFS)4&ii Since details of the method have been pub- 
lished elsewhere j 16 ' 17 ' 18 we now briefly discuss only the main 
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steps of the method. We first construct the limited func- 
tional space by starting from a Neel state with one hole with 
a given momentum k and zero phonon degrees of freedom 
l^ 0,0 ') = Ck|Neel;0), and applying the generator of states 

{l4f" M) >} = (tfki„ + fff r^L ' '}, where H kin and H g 
represent the first and the third term respectively of Eq. Q] 
This procedure generates exponentially growing basis space 
of states, consisting of different shapes of strings in the vicin- 
ity of the hole with maximum lengths given by Nh as well as 
phonon quanta that are as well located in the vicinity of the 
hole, at a maximal distance Nh- Parameter M provides gen- 
eration of additional phonon quanta leading to a maximum 
number N^ ax = MNh- Full Hamiltonian given by Eq. Q]is 
diagonalized within this limited functional space taking into 
account the translational symmetry while the continued frac- 
tion expansion is used to obtain dynamical properties of the 
model. The method treats spin, charge as well as lattice de- 
grees of freedom on equal footing. 
We define OC per doped hole 19 

*{u>) = ^ «r) - X ( W )) (2) 

POO 

XH=*/ e iw+t ([j(t)J(0)])dt (3) 
Jo 

where t = J2(ij) s s Ci, s represents the stress 

tensor, j = s s^m is the current operator, 

tij = —t for next nearest neighbors only and zero otherwise, 
and Rij = Rj — R^. We also note that in the case of next- 
neighbor tight binding models, (r) is related to the kinetic 
energy, (t^) = -(H kin )/2. 



m. CHARGE STIFFNESS AND SUM-RULES 

Charge stiffness per doped hole can be on a square lat- 
tice for the f-J-Holstein model computed via its spectral 
representation- 

= -\{0\H^\0) lE (0| ;f^ |0> , (4) 

= S tot -S^, (5) 

where S tot represents normalized optical sum-rule 
IZo a lA^ duj = 27r5 " ot ' while sre9 is defined by 
J Q+ a' „(a;)dw = irSfs and c'(cj) represents the real 
part of the optical conductivity tensor in Eq|2l We have 
computed in the single-hole ground-state, i.e. at 

k = (±7r/2, ±7r/2). It is well know that the dispersion E(k) 
is highly anisotropic around its single-hole minimum, which 
is in turn reflected in the anisotropy of the effective mass 
tensor 1 ' 8,21 . It is thus instructive to compute tensors repre- 
senting the charge stiffness as well as the OC in the direction 
of their eigen-axis, i.e. along the nodal ({tt/2, ir/2) — > (0, 0)) 
direction that gives (w), and along the anti-nodal 

((n/2, tt/2) — > (7T,0)) direction that leads to D±, and ct±(uj). 



In Fig. [If a) we present the charge stiffness vs. J/t for vari- 
ous values of EP coupling strength. To obtain accurate results 
in the strong EP coupling (SC) limit, we had to rely on only 
Ngt = 9786 different combinations of spin-flip states, while 
the total number of states, including phonon degrees of free- 
dom, was Ngt — 9 x 10 6 . To test the quality of A = results, 
we show with the dashed line Dm computed with zero phonon 
degrees of freedom using A^ st = 5 x 10 6 . Agreement with 
the A = case, obtained with N st = 9786 is rather surpris- 
ing, given the fact that results were computed using Hilbert 
spaces that differ nearly three orders of magnitude. This fast 
convergence is in contrast to calculations on finite-size clus- 
ters where due to the existence of persistent currents D varies 
rather uncontrollably between different system sizes 22 . 




I ■ i ■ i ■ i ■ i ■ i ■ i ■ i ■ I 

0.2 0.4 0.6 0.8 1 1.2 1.4 j/ t 



Figure 1: (Color online) a) Charge stiffness Dm ( D± in the inset), 

b) optical sum-rule S tot ({iV p h) = X^i a l a i m me insert) vs. J/t 
at u>o/t = 0.1 (full lines) and cuo/t = 0.3 (dashed lines). In this 
and all subsequent figures (except for the dotted line in a) or else 
otherwise indicated) we used: N h = 8, M = 7, and N st ~ 9 x 10 6 . 
Dotted line in a) was for A = obtained with states with zero phonon 
degrees of freedom and the following set of parameters: Nh = 14, 
M — 0, and Net ~5x 10 6 . 

Exploring further A = results we observe Dm ~ at 
J/t ~ indicating strong scattering on spin degrees of free- 
dom. With increasing J/t Dm steeply increases and around 
J/t ~ 1 reaches at Dm ~ 0.3 a broad maximum that as well 
coincides with the maximum of the bandwidth W— . In con- 
trast, the optical sum rule S tot — —{H^^/A monotonically 
decreases in the range < J/t < 1.0 with increasing J/t, 
Fig. |TJb). Due to strong anisotropy in E(k), D± remains 
nearly an order of magnitude smaller than Dm for J/t > 0.2 
(see the inset of Fig. [TJ a)). 

Turning to finite A, Dm expectedly decreases, due to addi- 
tional scattering on lattice degrees of freedom. The effect of 
A on the value of Dn however varies with J/t. This is best 
seen in the case of A = 0.25 and uJo/t = 0.1 where Dn is 
approximately equal to its A = value for J/t < 0.1, it then 
decreases with increasing J/t, reaching its minimum value 
around J/t ~ 0.5 and finally, for larger values of J/t > 0.8, 
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steeply increases. This non-monotonous begavior is as well 
reflected in the bell shaped average phonon number (N p h) vs. 
J/t, presented in the insert of Fig.[TJb). This behavior is also 
consistent with the non-monotonous functional dependence of 
X c (J/t), representing the crossover EP coupling strength to 
the SC regime, Refs.-^. 

We now make some general comments about the effect of 
the EP interaction on the correlated system at the onset of the 
SC regime. At small values of J/t < 0.1 EP coupling is less 
effective, which seems to be in contrast to naive expectations. 
We attribute this disentanglement from lattice degrees of free- 
dom to the increase of the kinetic energy and the vicinity of 
the Nagaoka regime. This effect particularly evident from 
the J/f-dependence of the average phonon number (-/V p h) at 
A = 0.3 (see the insert of Fig. QIb)) where an increase fol- 
lowed by a sharp drop of (iV p h) is seen with lowering of J/t. 
At the onset of the SC regime, i.e. at A ~ 0.25, EP coupling 
is most effective in the physically relevant J/t ~ 0.3 — 0.4 
regime, where there is a strong competition between kinetic 
energy and magnetic excitations. The critical A c as well 
reaches its minimum around J/t ~ 0.3 as shown in Ref. 16 . 
At larger J/t ~ 1 EP coupling becomes again less effective 
due to more coherent quasiparticle motion as reflected in the 
enhanced charge stiffess, quasiparticle weight, as well as the 
bandwidth 1 ^. 

The optical sum-rule S tot , presented in Fig. \V[b), as well 
decreases with increasing A. It however remains finite even 
deep in the SC regime where Dm ~ since 5 tot includes 
both coherent as well as incoherent transport. The latter re- 
mains finite due to processes, where the hole hops back and 
forth between neighboring sites while leaving lattice deforma- 
tion unchanged. Despite charge localization we thus expect 
nonzero optical response <r(u>) even deep in the SC regime, 
with its spectral weight shifted towards larger li and zero con- 
tribution at uj = 0. Due to localization we also expect OC to 
be isotropic in the SC regime, i.e. ctm (uj) ~ cr±(u>). 
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Figure 2: (Color online) a) Expectation value of the square of the 
electrical current along the nodal direction (j^) (S 1 ^ in the inset), 
b) expectation value of the square of the electrical current along the 
anti-nodal direction (5"j cg in the inset) vs. J/t at u>o/t — 0.1 
(full lines) and u)o/t = 0.3 (dashed lines). 



In the insets of Fig. [2] we show S™ & and S± s representing 
the integrated regular part of the OC. S™ g steeply decreases 
with increasing J/t due to the simultaneous increase of co- 
herent transport, captured by Dm, as well as due to decrease 
of S tot , see also Eq. [5] We observe, that EP coupling have 
little effect on S^° s for J/t < 0.4, since its value is rather in- 
dependent on A except deep in the SC regime, i.e. at A = 0.3 
in this particular case. Due to small values of D± we find 
greg _ 5 tot ) see Figs rjjand UJ. 

In Fig. |2] we present the average of the square of 
the electrical current defining the following sum-rule 
uia'^ „(u/)dcj = „), that furthermore represents the 
fluctuation of the current operator. In the ground state there is 
no persistent currents that usually appear on finite-size clus- 
ters, since our method is defined on an infinite lattice. This 
enables more reliable calculation of the charge stiffness. At 
A = (jf) (in Fig. Ha)) and (f x ) (in Fig. Hb)) display 

rather distinctive J/t dependence. While shows weak 

non-monotonous dependence on J/t, shows a substan- 
tial increase. With increasing A current fluctuations as well 
increase in both directions even though the increase is more 
pronounced in the case of In the SC regime we obtain 

(jf\) ~ Ox) as a consequence of localization due to lattice 
degrees of freedom. 
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Figure 3: (Color online) cry in a), b), and c), and u± in d), e), and f) 
for three different values of J/t as indicated in the figures for the t-J 
model for a single doped hole at A = and k = (n/2, 7r/2). Hilbert 
space with no phonon degrees of freedom and N Bt = 5 x 10 6 was 
used in all cases except in b) and e) where for comparison we in 
addition present calculations with N s t = 9786. In a), b), and d) 
we also show o xx = (ctm + u±)/2 using turquoise (dark grey) fill. 
Arrows in a) and b) indicate positions of lowest-energy peaks, ui\. 
Dashed lines in c) and f) are given by <r M ^ (ui) = tt/(zu>). Insert in c) 
represents scaling of Ui vs. J/t. In this and in the subsequent figures, 
the Drude peak is not shown. Artificial broadening e = Q.lt was 
used. Dashed areas in c) and f) delineate small frequency regimes 
(cu/t < 0.2) where at J/t = 0.05 EDLFS does not lead accurate 
results due to the vicinity of the Nagoka regime. 



IV. DYNAMIC PROPERTIES 

Turning to dynamic properties we first establish numerical 
efficiency of our method by presenting optical properties of 
the t-J model. In Fig. [3] we display different components of 
the conductivity tensor cr^^ui) in the single-hole minimum 
k = (7r/2, 7r/2), computed using EDLFS. At physically rel- 
evant value J/t = 0.4 we reproduce well known features, 
characteristic of a xx (uj): a) in the regime 1.6J < uj < 2t 
we find peaks forming a rather broad band, appearing within 
the well known mid-infrared (MIR) frequency regime, sepa- 
rated from the Drude peak (not shown) by a gap of the order 
of J and b) there is a broad featureless tail, extending to large 
frequencies, u> > It. MIR peaks for J/t > 0.2 scale with 
the exchange coupling (J/t) n , where r\ ~ 1. We stress that 
such scaling is consistent with local magnetic excitations as 
well as spin waves. Obtained scaling is however not consis- 
tent with the string picture where r\ = 2/3 (see also Refill). 
At J/t = 0.4 as well as at J/t = 1, the lowest peak appears 
at J/t ~ 1.6. Location of the lowest-frequency peak (indi- 
cated by arrows in Figs [3|a) and (b)) is surprisingly close to 
the location of the peak in OC of the t-J z model in the limit 
J z /t —>■ 00, given by o(uj) ~ t 2 /J z S(tu — §J Z )^. In this 
trivial case the peak appears at the frequency that corresponds 
to the energy (measured from the Neel state) of a single spin- 
flip, attached to the hole, created as the hole hops one lattice 
site from its origin in the undisturbed Neel background. It is 
somewhat surprising that such a naive interpretation seems to 
survive even in the (spin) isotropic t-J model and at rather 
small value of J/t = 0.4. The scaling of the position of 
the low-frequency peak closely follows the following expres- 
sion uji = 1.62( J/i) 108 , indicated by a dashed line, connect- 
ing the circles shown the insert of Fig. EJc). Our results of 
OC qualitatively agree with those, obtained on small lattice 
systems^. 

When conductivity tensor cr(uo) is computed in its eigen di- 
rections, distinct (incoherent) finite-w peaks are obtained in 
the case of cr\\(u>) and (t±(lu), as best seen at J/t = 1 and 
J/t = 0.4. For comparison we present in Figa. |3](a),(b), and 
(d) a xx (uj), that consists of all the peaks characteristic for both 
<j|| (uj) as well as a±(uj). The reason is, that the ground state 
at k = (tt/2, 7r/2) or S- point belongs to an irreducible repre- 
sentation Ei of the small group of k, i.e. C2. Current opera- 
tors j|| and j±, defining <r^ (uj) and <t±(uj) through EqsfJ] and 
[3]transform as distinct irreducible representations Ei and E2. 
Selection rules allow only transitions into states that transform 
according to a direct product of irreducible representations of 
the group C2. Since j x does not transform according to irre- 
ducible representations of C2, the above mentioned selection 
rules do not apply. 

In Figs. |3](b) and (e) we present as well results, computed 
on a much smaller set of states, i.e. with N st ~ 9786. Apart 
for a small shift of one of the MIR peaks at larger uj, the 
agreement with results, obtained with more than three orders 
of magnitude larger systems (N st — 5 x 10 6 ) underlines the 
efficiency of our method. Obtaining relevant results for the 
pure t-J model at moderate number of states is of crucial im- 
portance for successful implementation of additional lattice 
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Figure 4: (Color online) a xx for u)o/t — 0.3 in a), ujo/t = 0.1 
in b and c) at J/t = 0.3, and at k = (h/2,it /2). Total number 
of functions was Nat = 9 x 10 6 . Up to 56 phonon quanta was 
used to obtain accurate results for A > 0.2. Arrows in b) indicate 
uju = 16At where t = 0.45t. Units of a are arbitrary, yet chosen 
identical in a) and b); a different scale was used for c), nevertheless 
identical among different plots in c). Artificial broadening was set to 
e/t = 0.05. 



degrees of freedom. Last, we present in Figs. [3jc) and (f) 
results at small J/t = 0.05. Dashed lines represent known 
analytical estimate <j(uj) = it/(zuj), where z — 4, Ref. 26 . 
This result is characteristic for systems with a nearly constant 
density of states and diffusive hole motion where current ma- 
trix elements |(0|j M |n)| are roughly independent of n 26 . Good 
agreement with the analytical results in the small J/t limit 
is of particular importance since our method is by construc- 
tion, based on the existence of the long-range Neel order, tar- 
geted to be valid predominantly in the regime of intermediate 
to large values of the exchange constant J/t. We also note 
that in the limit of small- J/t optical properties for uj/t > J/t 
become isotropic. 

We now focus on the influence of increasing EP cou- 
pling A on optical properties of the t-J Holstein model in 
the adiabatic regime, i.e. for ujo/t = 0.3, Fig. |4ja) and 
uJo/t = 0.1, Fig.|4|b) and (c) the latter value being relevant 
for cuprates. Increasing EP coupling A leads to three main 
effects: a) the spectra progressively shift towards higher fre- 
quencies while the total spectral weight decreases (see also the 
inset of Fig.[3]T))), b) magnetic excitations that form a band in 
the MIR regime broaden and diminish with increasing A; they 
finally disappear in the SC regime where they are replaced 
by a broad polaron-like band that clearly originated from the 
renormalized MIR peaks. The peak of the well formed broad 
band at uJo/t = 0.1 in the regime 0.25 < A < 0.4 roughly 
scales with uj\\ ~ 16At where t — 0.45i represents renormal- 
ized hopping due to EP interaction. At uJo/t = 0.3 a broader, 
featureless band is formed, and c) a large gap opens in the SC 
regime. 

In contrast to numerical results of Ref. 13 , we observe no 
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Figure 5: (Color online) a M at A = 0.3; a) uja/t = 0.1, b) 
uo/t = 0.3, and A = 0.0 in c) at k = (ty/2,ty/2). Arrows in b) 
and c) indicate positions of the lowest-frequency state as it appears 
at respective values of J/t at A = 0. Artificial broadening was set to 
e/t = 0.05. 



separate peak at or slightly above the phonon frequency. This 
is more clearly seen in Fig. 0c) where (t xx (uj) is shown 
in an expanded frequency range. This result is consistent 
with DMFT calculations of Ref. 15 . Nevertheless, we find 
quantitative agreement at A ~ 0.24 with measurements on 
(Eui_ K Caj;)Ba2Cu306 in the low hole-doping regime pub- 
lished in Refrii In our calculation A ~ 0.24 represents the 
maximum EP coupling constant where the low-w peak, lo- 
cated at u>i ~ 1.56 J ~ l87meV ( choosing t = 400meV 
and J/t = 0.3), is just barely visible. This peak is, as dis- 
cussed above, due to the local magnetic excitation and re- 
mains separated from the continuum forming the rest of the 
MIR band. Experimental value of the corresponding peak is 
Wj Xp = l74meV—. The higher u— peak at u>n ~ 1.4t = 
560meV (experimental value is lu^ p — 590meV) corre- 
sponds to MIR band, slightly broadened and renormalized by 
phonon excitations. This part of OC is in agreement with cal- 
culations in Reffi^. 

Our explanation of the experimental results relies on the 
conjecture that lightly doped (Eui_ K Ca a; )Ba2Cu306 com- 
pound lies in the crossover from from weak to strong coupling 
electron-phonon regime where physical properties (quasipar- 
ticle weight, charge stiffness and dynamic properties) are ex- 
tremely sensitive to small changes of A. This is evident from 
Fig- E]and from results, published in Ref. 16 . MIR peak in OC 
is at A = centered around u>u = 2 J = 240meV. This value 
corresponds to the peak of the magnon density of states^ it 
however underestimates the position of the main peak, seen in 
the experiment of Ref. 13 . Increasing A beyond the weak cou- 
pling regime A > A c , the center of MIR peaks starts moving 
towards higher frequencies and broadens as it transforms into 
a wide polaron band, thus approaching the experimental value. 
Simultaneously the peak due to the local magnetic excitation 



at wj as well broadens and disappears above A > 0.24. 

The lack of a peak at ui ~ uiq in OC can be explained in 
simple terms in the large- J z /t limit of the simplified t-J z - 
Holstein model. Starting from a hole in the Neel background, 
the lowest energy contribution to a xx {uj) comes from the hop 
of the hole to the neighboring site. This move generates a sin- 
gle spin-flip with the energy E\ = 3J 2 /2 above the ground 
state. The contribution to OC that would include a single 
phonon excitation would thus be located at uj > 3 J 2 /2 + u)q. 

In order to explore the interplay of magnetic and polaronic 
degrees of freedom in the structure of <j xx (us) in more detail, 
we present in Figs. [2 a) and (b) comparison of optical spectra 
at fixed A = 0.3 and different values of the exchange inter- 
action J / 1. Decreasing J / 1 leads to a shift of the broad pola- 
ronic peak towards smaller values of u>. At smaller ojo/t = 0.1 
more pronounced structure abruptly appears at low ui/t < 0.5 
at small J/t = 0.05, Fig. 0a). At larger value of u Q /t = 0.3, 
Fig- 0b), a shoulder starts appearing at J/t = 0.3 in the low- 
uj regime that corresponds to the onset of the respective mag- 
netic peaks (as indicated by arrows in Fig. 0c)) of the pure 
t-J model. Below J/t < 0.2 well formed peaks emerge be- 
ing clearly of the magnetic origin. The disentanglement of 
lattice degrees of freedom, clearly seen in Fig. 0b), is consis- 
tent with DMFT calcutions 15 . 



V. SUMMARY 

In summary, we have explored effects of magnetic as well 
as lattice degrees o freedom on optical properties of the t-J- 
Holstein model. EDLFS captures well optical properties of a 
single hole in the t- J-Holstein model in the range of physi- 
cally relevant parameters of the model since it treats spin and 
lattice degrees of freedom on equal footing. Competition be- 
tween kinetic energy and spin degrees of freedom strongly in- 
fluences the coherent hole motion as measured by charge stiff- 
ness near the crossover to SC polaron regime. In the adiabatic 
regime increasing EP coupling leads to the shift of the OC 
spectra towards higher frequencies and broadening of peaks 
that in the pure t-J model originate in magnetic excitations. 
As an important as well as unusual finding we report a lack of 
a peak in the OC spectra at or slightly above the phonon fre- 
quency that we attribute to the inherently strong correlations 
that are present in the t-J model. This finding suggests that 
the two peak structure seen in recent optical measurements is 
entirely due to magnetic excitations. Based on our calcula- 
tions, the two peak structure can be explained with the obser- 
vation of local magnetic excitations, created by the hole mo- 
tion at lower frequencies and the contribution of spin waves, 
coupled via doped hole to lattice degrees of freedom at higher 
frequencies. 
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